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The electronic structure of Pbi_2,Sn2,Te is studied by using the relativistic Korringa-Kohn- 
Rostoker Green function method in the framework of density functional theory. For all concen¬ 
trations X, Pbi_2,Sn2,Te is a direct semiconductor with a narrow band gap. In contrast to pure lead 
telluride, tin telluride shows an inverted band characteristic close to the Fermi energy. It will be 
shown that this particular property can be tuned, first, by alloying PbTe and SnTe and, second, by 
applying hydrostatic pressure or uniaxial strain. Furthermore, the magnitude of strain needed to 
switch between the regular and inverted band gap can be tuned by the alloy composition. Thus, 
there is range of potential usage of Pbi_a:Sna;Te for spintronic applications. 


I. INTRODUCTION 

Among the narrow band gap semiconductors, SnTe 
and PbTe are in focus since long. The interest has been 
driven by their applications in electronic industry, essen¬ 
tially because of the effective control of charge carriers 
achieved in these systems. Tin and lead telluride alloys 
are attractive due to their use as infrared photodetec- 
tors^’^. Furthermore, PbTe and SnTe play an important 
role in the fabrication of thermoelectric materials^. Es¬ 
pecially by alloying, e.g. with bismuth (Bi) or antimony 
(Sb) a high figure of merit can be reached^’^. Recently, 
there has been revived interest in the system due to car¬ 
rier controlled ferromagnetism observed by suitable dop- 
ing6-8^ SnTe and its alloy with PbTe is found to be a 
topological crystalline insulator, where robust conduct¬ 
ing states are observed in their electronic structures®’ 

A topological surface state is also predicted for PbTe- 
SnTe multilayers^^. 

The topological characteristic in the electronic struc¬ 
ture is related to a fundamental band inversion at the 
Fermi level observed for tellurides like HgTe and SnTe. 
In the pure solids, both PbTe and SnTe have a direct 
band gap at the L-point of the Brillouin zone. However, 
whereas the analysis of the relativistic wave functions 
for PbTe show an even symmetry (L^) at the valence 
band edge and an odd symmetry {Lff) at the conduction 
band edge, the associated irreducible representations of 
the states at the band edges for SnTe are arranged in re¬ 
verse manner^®. As a result, the band gap of PbTe is said 
to be positive, whereas the band gap of SnTe is treated 
as negative. 

There exist various experimental reports on the change 
of the band gap with concentration x in Pbi_a;Sna;Te^®’^^. 
It is found that by increasing the concentration of Sn the 
band gap decreases linearly up to a certain concentration 
xq where the band gap vanishes. By further increasing 


X the band gap starts to increase linearly until the band 
gap of pure SnTe is obtained. In the past, the physical 
properties of the alloy have been studied theoretically by 
means of the k ■ tt perturbation method, where the band 
gap is treated as a function of the concentration x and 
the temperature^®’^®. Alloy properties were also studied 
by means of the tight-binding method^^’^®. An ab initio 
investigation of Pbi_a;Sna;Te was reported by Gao and 
Daw^®, where a quasirandom distribution of components 
within a supercell was used to study the alloy®®’®^. 

A conventional supercell approach is restricted by the 
number of atoms and only allows for particular con¬ 
centrations. To study the electronic properties of the 
Pbi_xSnxTe alloy within the whole concentration range 
0 > X > I the coherent potential approximation (GPA) 
was applied which has been used widely for the investi¬ 
gation of chemical and spin disordered alloys®®^®^. The 
main focus of this study is to investigate the size as well 
as the character of the band gap in dependence of the ap¬ 
plication of hydrostatic pressure and uniaxial strain. In 
contrast to hydrostatic pressure, the application of uni¬ 
axial strain is technologically more relevant in thin-film 
and multi-layers samples. 

The outline of the paper is arranged as follows: After 
introducing the computational methods in section II we 
discuss the electronic structure of pure PbTe and SnTe in 
section III A. The electronic structure of the Pbi_a;Sna;Te 
alloy as well as the influence of hydrostatic pressure is 
studied in section IIIB. Finally, the influence of uniaxial 
strain is discussed in section IIIC. 


II. COMPUTATIONAL METHODS 

The electronic structure calculations were performed 
by means of the Korringa-Kohn-Rostoker Green func¬ 
tion method (KKR)®"^ in the framework of density func- 
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tional theory^®’^®, as implemented within the computer 
code Hutsepot^’^. The atomic potentials of the scatter¬ 
ing centres were obtained self-consistently by applying 
the scalar relativistic approximation^®’^^ and the atomic 
sphere approximation^^. The evaluation of the Green 
function in terms of spherical harmonics was expanded 
up to a maximal angular momentum of / = 3. The en¬ 
ergy contour along the complex energy plane consisted 
of 24 points. As an approximation for the exchange- 
correlation functional we used the LibXC^^ implemen¬ 
tation of PBE^^, which is a generalized gradient ap¬ 
proximation. To properly simulate the Pbi_a;Sna;Te al¬ 
loy the coherent-potential approximation^^^^^ was ap¬ 
plied. Band structures were calculated in terms of the 
Bloch spectral function by using a fully-relativistic KKR 
method based on the solution of the Dirac equation. The 
relativistic KKR code was developed in order to describe 
the relativistic effects like spin-orbit coupling in a suit¬ 
able manner. Details can be found in Ref. 36. 

Pbi-a^Sn^^Te crystallizes in rock salt crystal structure. 
For the investigation of hydrostatic pressure on the elec¬ 
tronic structure a face-centred cubic (fee) primitive cell is 
chosen, consisting of one Te atom as well as one Pb or Sn 
atom, respectively. To improve convergence properties of 
the Green function, two additional empty spheres were 
added at non-occupied crystallographic sites. To simu¬ 
late uniaxial strain we chose a simple cubic (sc) unit cell, 
which contains four times as many atoms as the primi¬ 
tive fee cell. The lattice vectors are given by a = (a, 0,0), 
b = (0, b, 0) and c = (0, 0, c). The relation uq = bo = cq 
holds for the equilibrium lattice constants. The associ¬ 
ated dimensions of the fc-mesh within the Brillouin zone 
were 20 x 20 x 20 for the fee cell, as well as 12 x 12 x 12 
for the sc cell. 

The atomic sphere approximation is a coarse approxi¬ 
mation for estimating equilibrium lattice constants. We 
thus employed the Vienna ab initio simulation package 
(VASP) for estimating the lattice parameters of the pure 
systems PbTe and SnTe. For the calculations 8x8x8 
/c-point mesh was chosen. Galculations were done using 
the high performance mode using a plane wave energy 
cut-off of 350 eV. The lattice parameters of the alloy 
are calculated from Vegard’s law^^ which states that the 
lattice constant changes linearly for 0 < a; < 1. Test cal¬ 
culations with VASP using a large unit cell agreed to the 
approximation fairly well. 


III. RESULTS AND DISCUSSION 


A. Electronic structure of PbTe and SnTe 

The lattice parameters and the band gap for PbTe 
and SnTe obtained with VASP are tabulated in Tab. I. 
The equilibrium lattice constants are found to be slightly 
larger than those reported in experiments^®’^®. It was 
observed that spin-orbit interaction (SOI) does not in¬ 


fluence the lattice constants. The band gap of PbTe and 
SnTe are estimated to be 0.81 eV and 0.05 eV in GGA. 
With the incorporation of SOI the valence bands and 
conduction bands split. The splitted bands effectively 
move towards each other thereby reducing the band gap 
to 0.07 eV and —0.11 eV for PbTe and SnTe, respectively 
(refer Tab. I). Here, the negative band gap signihes the 
inversion of symmetry of the states at the band edges. 
Although the splitting for SnTe is smaller than for PbTe, 
the smaller non-relativistic band gap of SnTe allows band 
inversion. In agreement with previous theoretical stud- 
ies^®^'*® it was found that SOI is crucial for reducing the 
band gap for both compounds with characteristic band 
inversion for SnTe. 

In the framework of the fully relativistic Green func¬ 
tion method, the Bloch spectral function (BSF) A^^k, E) 
is calculated which is given by^® 


A-Q^k, E) — 


-Im Tr 

TT 


j d^r,air„r,+R„E) 

RjGC 


■ ( 1 ) 


The summation on the right-hand side of the equation 
extends over all vectors Rj of the lattice C. The trace 
is a summation of the diagonal components of the 4x4 
dimensional multiple-scattering Green function"^^ G and 
E denotes the energy. 

In the band structure of a regular semiconductor, the 
valence band is made up of anion character while the 
conduction band is made up of cation character. In band- 
inverted semiconductors this property is reversed. Thus, 
the difference of the site-resolved BSF from the cation 
and the anion. 


AAb = (2) 

can be used as a tool for analysis of band inversion prop¬ 
erties. In Fig. I a) and Fig. 1 b), AAb is plotted for PbTe 
and SnTe, respectively. From the intermixed colours near 
the Fermi level which is set to zero the band inversion 
characteristic of SnTe is clearly revealed. 


Table I. Equilibrium lattice parameters and band gap calcu¬ 
lated using VASP in comparison with experimental values. 
Calculations were performed with and without spin-orbit in¬ 
teraction, denoted as GGA and GGA+SOI, respectively. 


Method 


PbTe 


SnTe 

a (A) 

E, (eV) 

a (A) 

E, (eV) 

GGA 

6.565 

0.81 

6.408 

0.05 

GGA -b SOI 

6.562 

0.07 

6.410 

-0.11 

Expt. 

6.462®® 0.19^® 

0.32^® 

6.303®® 

-0.2—0.3'‘^ 
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Figure 1. The difference between the cation and the anion 
contribution within the Bloch spectral function (AAb) plot¬ 
ted for the T — L — W high-symmetric line in the Brillouin 
zone for PbTe and SnTe, respectively. The Fermi energy is 
set to zero. The red color signifies that the main contribution 
to the BSF is from the anion atom (Te), while the blue color 
signifies that the main contribution is from the cation atom 
(Pb or Sn). 


B. Electronic structure of Pbi_a;Sna;Te and 
hydrostatic pressure 

In our studies of the alloy we stick to the rock salt crys¬ 
tal structure over the whole range of concentration, even 
though a few reports on a temperature induced rhomb- 
hohedral phase exist^®“®°. The electronic structure of the 
Pbi_a;Sna;Te was investigated for various concentrations 
X and various lattice constants. The band gap energy 
was estimated by calculating the BSF [Eq. (1)] at the L- 
point of the Brillouin zone within a small energy region 
[Ep — A, Ep + A] (Aw 0.2 eV) about the Fermi energy 
Ep using Np = 700 energy points. All together the cal¬ 
culations were performed for 11 different concentrations 
and 20 different lattice parameters each, to produce a rea¬ 
sonable array of data. Since it is not possible to resolve 
the band gap in the region where the band gap changes 
sign, this array was interpolated using Mathematical^ to 
obtain a well defined curve. 

The result in terms of band gap versus concentration 
and lattice constant is illustrated in Fig. 2. The vertical 
dashed lines show the equilibrium lattice parameter of 
SnTe and PbTe. Following the line of Vegard’s law (an¬ 
gular dashed line) from pure PbTe to pure SnTe, it can 
be seen, that with increasing amount of Sn the band gap 
energy decreases until it vanishes for a Sn concentration 
of « 40% (vanishing intensity of red shade). This is in 
good agreement with experimental reports^^’®^’®^. Low¬ 
ering the Pb concentration further increases the band 
gap, but with an inverted band characteristic at the L- 
point (growing intensity of blue shade). Hence, the size 
of the band gap as well as the band characteristic of 
Pbi_a;Sn 2 :Te can be tuned with changing the concentra¬ 
tion of Pb. 

Starting again with pure PbTe and leaving the con¬ 
centration constant, a similar behaviour can be seen for 
compressing the crystal. Again, compression leads to a 
decrease of the band gap until the band gap vanishes at a 


lattice constant of about 6.43 A. This value corresponds 
to a lowering of the lattice constant of about 2%. By 
using the Murnaghan equation of states, 


P{V) 



( 3 ) 


together with the values Bq = 37.73 GPa and B' = 3.52 
obtained from our calculations, which are in good agree¬ 
ment with reported values for PbTe®'^, it can be verified 
that these conditions can be obtained by applying a pres¬ 
sure of Pc « 2.63 GPa. Further compression will open 
the band gap again but with inverted band characteris¬ 
tic. Our estimation of the critical pressure Pc is in very 
good agreement with the experimental value of 3 GPa®^. 
According to Fig. 2, Pc can be tuned by varying x in the 
Pbi_ 2 ;Sn 2 .Te alloy. 

In order to understand the electronic effects for the 
band gap change for the Pbi_a;Sna;Te alloy, we refer to 
the expression of the band gap achieved by means of the 
tight-binding model®®, 

F;g«Ao + 104(Ar'-A2-i)/3, (4) 

where the abbreviations 


Aq ^p^Pb/Sn ^p,Te: 

^ = Cp.Pb/Sn — Cs.Te , and 

A 2 = ^p,Te ^s,Pb/Sn¬ 
are used. Here, e^x and ei^x are the orbital eigenvalues 
obtained by a non-relativistic and a relativistic treatment 



6.3 6.4 6.5 6.6 6.7 


Lattice constant (A) 

Figure 2. Band gap as a function of the lattice constant and 
the concentration x in Pbi-jjSnjjTe. The gradient shade with 
red color represents regular band gap and the blue color rep¬ 
resents the inverted band gap. The dashed vertical lines rep¬ 
resent the equilibrium lattice parameters of SnTe and PbTe 
and the straight dashed line across represents the lattice pa¬ 
rameters versus concentration according to Vegard’s law®^. 
The dot represents Eg = Q along the line obeying Vegard’s 
law. 
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Figure 3. Band gap of Pbi_i:Sni:Te as a function of c/co 
obtained from fully-relativistic Green function calculations. 


for the quantum number I of element X. They are related 
by, £; X = e;_xi2Ax, where Xx is the spin-orbit coupling 
strength of element X. In principle, there are two com¬ 
peting effects which influence the band gap, 1) change 
in Aq and t through the volume of the unit cell and 2) 
the change of the SOI-strength A through the concentra¬ 
tion X. It has been suggested that the former criterion is 
decisive for the band gap change versus concentration®®. 

From the results we conclude that both effects cause 
a linear change of the band gap energy. In fact, a lin¬ 
ear scaling of the band gap among the values —0.11 eV 
to 0.07 eV for X = 1 and x = 0 gives zero band gap 
at X = 0.39, which is very close to the value estimated 
from KKR CPA calculations, x = 0.4. Our conclusion 
is in agreement with Wexler and Morgan®^ who studied 
Pbi_a;Sn 2 :Te using a strain field and including relativis¬ 
tic effects in the framework of an augmented plane-wave 
method. It must be noted that the SOI modeled sepa¬ 
rately also leads to an almost linear change of the lattice 
constant®®. In accordance to Ref. 56, this emphasizes 
that the band inversion characteristics can be controlled 
by changing the volume. 

In summary, a linear behaviour for the change of the 
band gap can be observed along the line of Vergard’s 
law, where the equilibrium lattice constant changes lin¬ 
early with respect to the composition of the alloy. The 
same holds true for the change of the band gap for com¬ 
pressed cells below the equilibrium volume but parallel 
to Vegard’s line. However, for larger volumes (parallel to 
Vegard’s line), the SOI induces nonlinear effects which 
are reflected as the white region in Fig. 2, where Eg = 0 
appears nonlinearly. 


C. Strained Pbi_i,Sna;Te alloy 

Apart from the effect of hydrostatic pressure on solids, 
uniaxial and biaxial strain in the alloy is technologically 
more relevant. Thus, we have undertaken the study of 
band inversion characteristics in Pbi_ 3 ;Sn 2 ,Te upon the 
action of uniaxial strain. Our attention is towards the 


physically interesting concentration range of 60 % Pb and 
higher (refer to Fig. 2). Therefore, the concentrations 
100 %, 90 %, 78 % and 60 % were chosen to investigate 
how uniaxial strain influences the band gap. The unit cell 
was compressed or strained along the crystallographic c- 
direction. The value of a (= b) is thereby derived with 
the assumption that the volume of the unit cell is fixed as 
c is changed. This approximation is reasonable, since the 
change of the volume is sufficiently small with respect to 
the concentration x. 

The band gap as a function of c for the above men¬ 
tioned concentrations is shown in Fig. 3. It can be 
verified that both, compression as well as strain leads 
to a lowering of the band gap with a parabolic depen¬ 
dence of c/cq. In both cases, band inversion can be trig¬ 
gered if the change of the lattice constant c is reason¬ 
ably large. For pure PbTe the band gap vanishes for the 
values c/cq « 0.96 and c/cq « 1.045. By alloying, the 
band gap at the equilibrium changes {a = b = c = cq) 
and therefore the transition between inverted and non- 
inverted band characteristic already occurs for smaller 
deformations. For a Pb concentration of 60%, the un¬ 
strained Pba;Sni_ 3 ,Te is metallic. Our results show that 
by applying tensile or uniaxial strain it is possible to in¬ 
duce an inverted band gap into the alloy. 

This conclusion is in agreement with the study of 
Qian et al. who did a similar investigation for bulk 
SnTe®®. Furthermore, they proposed that strained SnTe 
membranes can be used as frequency resolved infrared 
photodetectors. By choosing samples of Pba;Sni_ 2 ,Te and 
gradually changing the composition a comparable device 
can be constructed which is less fragile in comparison to 
thin SnTe membranes. 

Two physical parameters, namely, strain and alloy 
composition allow tuning the band gap and the spin- 
texture of band edge states. Potential applications of 
this alloy lie in switching between a topological and a 
regular semiconductor by applying strain in spintronic 
devices. Furthermore, the effective strain can be con¬ 
trolled by choosing a proper composition of the alloy. 
Thus, besides the application as infrared photodetectors 
for a wide range of frequencies the tuning of band in¬ 
version leads to exciting possibilities for novel functional 
devices such as spintronic-switches. 


IV. CONCLUSIONS 

We have studied the band inversion characteristics of 
Pbi_ 2 ;Sna;Te using the relativistic KKR method com¬ 
bined with the CPA. In agreement with experiments it is 
found that an inversion of the band characteristic at the 
L-point occurs for a lead concentration of 60 %. By ap¬ 
plying hydrostatic pressure or uniaxial strain, the band 
gap can be tuned as well. Whereas the change of the 
band gap is more or less linear close to the equilibrium 
lattice constants for hydrostatic pressure, the change of 
the band gap shows a parabolic behaviour for the ap- 
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plication of uniaxial strain. Hence, for the latter, band 
inversion can be obtained by either compressing or strain¬ 
ing the sample along the crystallographic c-direction. 

As a consequence, the application of strain or hydro¬ 
static pressure to Pbi_ 2 ,Sn 2 ;Te allows to switch between 
a topological and a regular semiconductor. In this con¬ 
nection, the magnitude of the needed external pressure 
can be controlled by changing the composition of the al¬ 
loy. Therefore, Pbi_a;Sna;Te is a prominent candidate for 
the fabrication of topological-switches. 
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